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1  Introduction 


Let  V  be  a  three-dimensional  real  inner  product  space,  and  let  R  denote  the  real 
numbers  or  scalars.  Let  Sym  denote  the  linear  space  of  symmetric  (second-order) 
tensors  on  V.  Let  U  he  a,  nonempty  subset  of  Sym  which  is  invariant  in  the  sense 
that  A  €  ZY  iff  QAQ^  G  U  for  each  orthogonal  tensor  Q.  Then  a  scalar-valued 
isotropic  function  on  ZY  is  a  function  ^  :  Z/  R  with  the  property  V’(QAQ^)  =  ■0(A) 
for  each  A  eU  and  each  orthogonal  tensor  Q.  A  tensor-valued  isotropic  function  on 
Id  is  a,  fvmction  ^  :  Id  Sym  with  the  property  $(QAQ^)  =  Q$(A)Q^  for  each 
A  Eld  and  each  orthogonal  tensor  Q.  For  any  symmetric  tensor  A,  let  *A  E  {1, 2, 3} 
denote  the  number  of  distinct  eigenvalues  of  A.  Let 

Sym„  :=  {A  E  Sym  :  *A  =  n}  (n  =  1,2, 3) .  (1.1) 

In  particular,  Sym^^  is  the  one-dimensional  subspace  of  all  spherical  tensors  al,  a  €  R. 
Each  set  Sym^  is  invariant,  and  Sym  =  Sym-^  U  Sym2  U  Sym^.  Let 

U„:=UC\Sym^  =  {AEU  ■.*A  =  n]  (n=:  1,2,3).  (1.2) 

Then  each  set  is  invariant,  and  ZZ  =  ZYi  U  ZZ2  UZZ3.  A  scalar- valued  or  tensor- valued 
function  on  Id  is  isotropic  iff  its  restriction  to  each  of  the  invariant  subsets  ZZ„  is 
isotropic. 

According  to  a  well-known  theorem  of  Rivlin  and  Ericksen  [1,  §29,  §39],  a  function 
$  :  ZZ  Sym  is  isotropic  iff  there  are  scalar- valued  isotropic  functions  a, on  Id 
such  that 

$(A)  =  a(A)I  -h  ^{A)A  -h  7(A)A2,  (1.3) 

where  I  denotes  the  identity  tensor.  Since  the  publication  of  Rivlin  and  Ericksen’s 
classic  paper  in  1955,  various  proofs  of  this  representation  theorem  have  appeared  in 
the  mechanics  literature;  cf.  Serrin  [2],  [3,  §59],  Truesdell  &  Noll  [4,  §12],  Wang  ic 
Truesdell  [5,  §111.2],  Gurtin  [6,  §37],  and  Wang  [7].^  The  proofs  in  [l]-[7]  have  the 
following  features  in  common.  The  authors  establish  (1.3)  by  showing  that  it  holds 
on  each  of  the  subsets  Id^^^ld^Ms-  particular,  they  show  that  a, y0,7  are  uniquely 
determined  by  $  on  ZZ3.  They  set 

7(A)  =  0,  VAsMj.  (1,4) 

and  show  that  there  are  unique  functions  a,  /?  on  ZZ2  satisfying  (1.3)  and  (1.4).  Finally, 
they  set 

/5(A)  =  7(A)=0.  VA€M..  (1.5) 

and  show  that  there  is  a  unique  function  a  on  ZZj  satisfying  (1.3)  and  (1.5).  Isotropy 
of  the  coefficients  is  established  in  one  of  two  ways:  either  by  verifying  indirectly  that 

^Applications  of  (1.3)  to  the  response  functions  of  isotropic  materials  are  discussed  in  [l]-[6]. 
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7(QAQ^)  =  7(A),  for  example,  without  actually  solving  for  the  coefficients  (cf.  [4]- 
[7]),  or  by  deriving  explicit  isotropic  formulas  for  the  coefficients.  The  latter  approach 
was  used  by  Rivlin  and  Ericksen  [1]  (cf.  also  Serrin  [2,  3]),  and  their  formulas  for  the 
case  #A  =  3  are  of  interest  here. 

Let  01,02,03  denote  the  eigenvalues  of  the  symmetric  tensor  A,  let 
denote  the  principal  invariants  of  A  ,  and  let  II a-,  HI  a  denote  the  second  and  third 
moments  of  A.  Then  (cf.  Ericksen  [8,  §38]) 


Ia  =  Oi  +  02 +  03  =  tr  A  =  A*I,  (1.6) 

IIa  =  O1O2  +  02O3  +  03O1  =  2(1^  ~  -^A.)  5  (l-^) 

IIIa  =  010203  =  det  A  =  \Il  -  ^Ia^Ia  +  k^A  ?  (1-8) 

^A  =  o^  +  ol  +  Og^  =  trA^  =  A-A  =  ||A|P  ,  (1.9) 

in  A  =  of  +  of  +  of  =  tr  A3  =  /^  -  37a/7a  +  3/74  •  (I-IO) 

Here  ||A||  denotes  the  norm  of  A  corresponding  to  the  inner  product  A  •  B  =  tr  (AB) 
on  Sym.  From  (1.3)  we  obtain  the  following  system  of  equations  for  q:(A),  /?(A),  7(A): 

[3  4  1  r  ^(A)  1  r  tr$(A) 

4  ^  777a  ^(A)  =  tr(#(A)A)  .  (1.11) 

_  74  IIIa  trA4  J  [  7(A)  J  [  tr($(A)A2)  _ 

Let  Ma  denote  the  coefficient  matrix  in  (1.11).  By  the  Cayley-Hamilton  theorem, 

A3  =  7a A2  -  77aA  +  IIIaI  •  (1.12) 


On  multiplying  (1.12)  by  A  and  taking  the  trace,  we  obtain  an  expression  for  tr  A^ 
which,  together  with  (1.6)-(1.10),  can  be  used  to  express  detMA  as  a  polynomial  in 
either  the  principal  moments,  the  principal  invariants,  or  the  eigenvalues  of  A.  In 
particular,  we  find  that  det  Ma  is  just  the  cubic  discriminant  Aa  of  the  characteristic 
polynomial  of  A: 

detMA  =  AA  =  187a77a777a-47^777a  +  7^2jj^- 477^ -27777^ 

=  (fli -02)3(02 -03)^(03-01)^  •  (1*13) 

Since  Aa  ^  0  iff  =  3,  (1.11)  has  a  unique  solution  for  q:(A),  ;5(A),  7(A)  iff  *A  =  3, 
in  which  case  Cramer’s  rule  yields  the  explicit  formulas  (Rivlin  and  Ericksen  [1,  §29]) 

ct(A)  =  det  mI'VAa  ,  ^(A)  =  det  M^/Aa  ,  7(A)  =  det  iWf /Aa  ,  (1.14) 

where  denotes  the  matrix  obtained  by  replacing  the  nth  column  of  Ma  by  the 
column  vector  on  the  right-hand  side  of  (1.11).  Then  the  isotropy  of  a,  ^,7  on 
follows  from  (1.14),  the  isotropy  of  the  principal  invariants  and  the  moments,  and  the 
isotropy  of  #. 


The  definitions  and  results  above  place  no  restrictions  on  U  other  than  invariance. 
However,  since  we  are  interested  in  the  relationship  between  the  smoothness^  of  $  and 
the  smoothness  of  the  coefl&cients  in  its  representations,  for  the  remainder  of  the  paper 
we  assume  that  the  domain  U  of  the  tensor- valued  isotropic  function  $  is  also  an  open 
subset  of  Sym.  Then  and  U—Ux  =  U  ^4  are  open  since  Sym^  and  Sym  —  Sym^  = 
Sym2  U  Sym^x  open.  And  since  Sym^  is  dense  in  Sym,  and  U  Uz  axe  dense  in 
U.  Now  from  (1.3)  it  follows  that  any  smoothness  properties  shared  by  a,/?, 7  on  an 
open  subset  U  of  U  are  also  shared  by  $.  Conversely,  since  the  principal  invariants 
and  the  moments  of  A  axe  smooth  functions  of  A,  the  formulas  (1.14)  imply  that  on 
the  open  subset  U3  the  coefficient  functions  a,  7  inherit  any  smoothness  properties 
of  $.  Thus  the  problem  of  determining  the  relationship  between  the  smoothness 
of  $  and  the  smoothness  of  the  coefficients  in  its  representation  (1.3)  reduces  to 
determining  the  smoothness  properties  that  a,  0, 7  inherit  from  #  at  points  in  the 
set  lA  —Uz  =  ZYi  U  Uz-,  which  is  closed  and  nowhere  dense  relative  to  TA^  Rivlin 
and  Ericksen  [1,  §29]  obtained  explicit  formulas  for  q;(A)  and  /?(A)  when  #A  =  2 
and  7  satisfies  (1.4);  and  on  ZYj  we  simply  have  a  =  |tr$  when  0  and  7  satisfy 
(1.5).  But  these  formulas  axe  of  no  use  for  the  problem  at  hand.  There  are  infinitely 
many  other  sets  of  isotropic  functions  a,/?, 7  for  which  the  representation  (1.3)  holds 
(though  all  such  sets  necessarily  coincide  on  Uf).  The  particular  set  of  coefficient 
functions  satisfying  (1.4)  and  (1.5),  while  convenient  for  the  proof  of  (1.3),  may  fail 
to  be  continuous  on  U  even  for  polynomial 

Now  suppose  $  is  continuous,  so  that  a,  ^,7  are  necessarily  continuous  on  Uz- 
Since  Uz  is  dense  in  U,  by  a  well-known  theorem®  of  analysis,  oc,0,j  have  unique 
continuous  extensions  to  all  of  U  iff  their  limits  from  within  Uz  exist  at  each  point 
of  UiD  Uz',  and  it  is  not  hard  to  show  that  when  these  continuous  extensions  exist 
they  are  isotropic  on  U  (cf.  Proposition  2.1).  However,  in  the  explicit  formulas  (1.14) 
for  a,0,^  on  Uz,  the  denominator  A^.  0  as  A  ^  U  Uz,  so  sufficient  conditions 

for  the  existence  of  these  limits  are  by  no  means  obvious.®  This  problem  was  first 

^For  the  present  discussion  we  do  not  attach  any  precise  meaning  to  the  term  “smooth”;  e.g.,  it 
could  mean  continuous,  differentiable,  or  ^-times  continuously  differentiable  (denoted  by  C*). 

®If  U  lies  in  the  set  of  nonsingular  symmetric  tensors,  then  on  multiplying  (1.12)  by  and 
using  the  representation  (1.3),  we  obtain  an  alternate  representation  which  is  widely  used  in  studies 
of  nonlinear  isotropic  elastic  solids  (cf.  [4]-[6]);  ^(A)  =  a(A)I ^(A)A  -t-  7(A) A~^,  where 

a(A)  =  a(A)  -  i7A,7(A),  0{A.)  =  /?(A)  -t-  4v7(A),  and  7(A)  =  IZiv7(A).  It  follows  that  a,  /?,  7  are 
smooth  at  a  point  A  iff  a,  jd,  7  are  smooth  at  A. 

^For  example,  by  (112)  we  see  that  the  isotropic  function  #(A)  =  A®  on  W  =  Sym  has  the 
representation  (1.3)  with  smooth  coefficient  functions  a(A)  =  EIx,  /?(A)  =  — JTa,  t(A)  =  /a, 
whereas  the  a,0,j  satisfying  (1.4)  and  (1.5)  are  discontinuous  since  they  must  agree  with  the 
preceding  solution  on  Sym^. 

®Cf.  (3.15.5)  in  Dieudonne  [9]. 

®It  is  easily  seen  that  continuity  of  #  is  not  sufficient.  Consider  the  continuous  isotropic  function 
$  on  Sym  defined  by  #(0)  =  0  and  #(A)  =  ||A||~^^^(A  -f  A^)  for  A  ^  0.  Then  0{A.)  =  7(A)  = 
l|A||  on  Sym^x  since  these  functions  are  unbounded  on  every  neighborhood  of  0,  they  cannot 
be  extended  continuously  to  Sym. 
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addressed  by  Serrin  [2],  who  showed  that  the  coefficients  a,/3,j  in  the  representaion 
(1.3)  can  be  extended  continuously  from  to  if  $  is  C^.  In  proving  the  existence 
of  the  limits  of  q:(A),  j3{A),^{A)  as  A  ->  B  €  U  Z/2  from  within  Z/j,  he  treated  the 
cases  B  €  2^2  and  B  €  Z/j  separately.  These  cases  are  referred  to  as  double  and  triple 
coalescence,  respectively.  For  double  coalescence,  his  proof  requires  only  that  $  be 
C^.  Serrin  [2]  also  gave  an  example  of  a  isotropic  function  $  on  Sym  for  which 
the  coefficients  cannot  be  extended  continuously  to  Sym-,  his  example  is  discussed  in 
Section  2.  There  appears  to  have  been  no  additional  contributions  to  this  problem 
until  the  paper  by  Man  [10]  thirty-five  years  later.  Man  remarks  that  Serrin’s  proof, 
particularly  for  the  case  of  triple  coalescence,  was  given  only  in  outline.  Man  gives  a 
different  proof  (with  full  details)  of  the  continuity  of  a,  7  for  this  case.  His  proof 
has  the  advantage  that  $  is  only  required  to  be  C^.  On  combining  these  results,  we 
have 

Theorem  1.1  (Serrin  [2]  and  Man  [10])  If  ^  isC^  onU  (respectively,  on 
IA2  U  ZY3J,  then  there  are  unique  continuous  scalar-valued  isotropic  functions  a,  7  on 
U  (resp.  U2  U  U^)  for  which  the  representation  (1.3)  holds. 

I  am  not  aware  of  analogous  results  for  the  differentiability  of  the  coefficients.’'  One 
of  the  goals  of  this  paper  is  to  fill  this  gap  in  the  literature. 

In  the  next  section  I  confirm  Theorem  1.1  using  a  method  which  is  substantially 
different  from  that  of  Serrin  and  Man.  Their  proofs  of  the  existence  of  limits  of  the 
coefficients  at  points  in  tii  U  IA2  are  complicated  by  the  fact  that  the  denominator 
in  the  formulas  for  these  coefficients  goes  to  zero.®  Instead,  in  Section  2  I  show  that 
on  Z/3  there  are  simple  formulas  for  the  coefficients  in  terms  of  the  first  or  second 
derivatives  of  From  these  formulas  it  follows  trivially  that  the  coefficients  have 

'’Of  course,  it  is  well-known  (cf.  [l]-[4])  that  if  $  is  polynomial  in  the  sense  that  the  components 
of  ^(A)  are  polynomials  in  the  components  of  A,  then  the  coefficients  in  (1.3)  may  be  expressed  as 
polynomials  in  the  principal  invariants  of  A;  cf.  Serrin  [2],  [3,  §60],  for  a  simple  proof.  Regarding 
polynomial  dependence  of  response  functions,  'ftuesdell  &  Noll  [4,  p.  61,  footnote]  remark  “To  us, 
assuming  polynomial  dependence  seems  not  only  unnecessary  ...but  unphysical.  We  see  no  sign 
that  nature  loves  a  polynomial,  and  polynomial  dependence  is  not  even  invariant  under  change  of 
strain  measure.”  The  primary  application  of  the  above  result  on  polynomial  dependence  would  seem 
to  be  to  the  approximation  of  tensor- valued  isotropic  functions.  If  #  is  C”’  then  its  rth-order  Taylor 
approximation  at  a  spherical  tensor  cl  is  an  isotropic  polynomial  of  degree  r  and  hence  has  a 
representation  of  the  form  ^(''^(A)  =  ar(A)I  -|-  y0r(A)A  -f  7r(A)A^,  where  ar(A),  /?r(A),  7r(A) 
are  polynomials  in  the  principal  invariants  of  A.  However,  this  does  not  imply  that  Pr,  7r 
rth-order  approximations  of  a,  /?,  7  at  cl.  Indeed,  example  (2.2)  in  the  next  section  shows  that  for 
the  case  r  =  1,  the  polynomial  functions  ai,  ^i,  71  may  be  identically  zero  even  if  a,  /?,  7  take  on 
all  real  values  in  every  neighborhood  of  every  spherical  tensor. 

Serrin  [3,  p.  235,  footnote  2]  states  that  the  theorem  on  polynomial  dependence  may  be  generalized 
to  analytic  $,  the  coefficients  in  this  case  being  analytic  functions  of  the  principal  invariants. 

®They  worked  with  explicit  formulas  for  o;,/?,  7  on  U3  similar  to  (114)  but  expressed  in  terms  of 
the  eigenvcilues  of  A  and  ^(A).  These  formulas  were  noted  by  Rivlin  and  Ericksen  [1,  §29];  cf.  also 
I^uesdell  &  Noll  [4,  §48]. 
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limits  at  points  in  ZYj  U  IA2  from  within  certain  subsets  of  and  then  existence  of  the 
limits  along  arbitaraxy  paths  in  Z/j  can  be  inferred  from  the  isotropy  of  the  coefficients 
(cf.  Proposition  2.1).  The  proof  does  not  require  a  separate  treatment  of  double  and 
triple  coalescence.  A  similar  procedure  is  used  in  Section  4,  where  I  establish 

Theorem  1.2  7/  $  is  on  U,  then  the  continuous  isotropic  coefficients 
in  the  representation  (1.3)  are  on  U.  If  $  is  (resp.  C^)  on  U2  U  U^,  then 
a:,/?, 7  are  (resp.  C^)  on  ZY2  U  7/3. 

My  approach  also  differs  from  Serrin’s  and  Man’s  in  that  I  utilize  the  decomposi¬ 
tion  of  #  into  its  spherical  and  deviatoric  parts.  For  any  symmetric  tensor  A,  let  A* 
denote  the  deviatoric  part  of  A: 

A*  :=  A  -  |7aA  .  (1.15) 

Then  A  is  deviatoric  if  A  =  A*  (equivalently,  Xv  =  0);  A  is  spherical  iff  A*  =  0. 
Let  Sym*  denote  the  subspace  of  Sym  consisiting  of  all  deviatoric  symmetric  tensors. 
The  deviatoric  part  of  $  is  the  function  :U  Sym*  defined  by  ^*(A)  :=  #(A)*. 
Then 

$(A)  =  i(tr$(A))I-F$*(A);  (1.16) 

and  since  tr  :  Sym  R  and  (  )*  :  Sym  — >  Sym*  are  linear,  the  isotropic  functions  tr  $ 
and  are  both  smooth  iff  $  is  smooth.  There  are  several  advantages  to  using  the 
decomposition  (1.16):  it  simplifies  some  of  the  computations;  it  allows  us  to  obtain 
the  conclusions  in  the  smoothness  theorems  under  slightly  weaker  conditions  on  $, 
e.g.,  the  conditions  on  #  in  Theorem  1.1  can  be  replaced  by  the  same  conditions  on 
together  with  the  continuity  of  tr  it  reveals  that  relatively  stronger  smoothness 
properties  hold  for  some  of  the  coefficients  in  other  useful  representations  for  $  and 
(see  below);  and  it  more  clearly  reveals  the  nature  of  the  possible  loss  in  smoothness 
of  the  coefficients. 

By  taking  the  deviatoric  part  and  the  trace  of  the  representation  (1.3),  we  obtain 

$*(A)  =  ^(A)A*  +  7(A)(A2)*  (1.17) 

and 

tr  $(A)  =  3a(A)  +  4/3(A)  +  .  (1.18) 

On  substituting  (1.17)  (or  (1.20)  below)  into  (1.16)  we  obtain  a  representation  for  $ 
which  is  similar  to  (1.3)  but  has  the  property  that  the  coefficient  of  I  is  as  smooth 
as  Although  are  not  uniquely  determined  by  from  (1-18)  we  see  that  a 

is  uniquely  determined  by  tr$,  /?,  and  7  and  that  a  has  any  smoothness  properties 
shared  by  tr  #,  and  7.  Therefore,  it  remains  to  establish  the  smoothness  of  /?  and 
7  at  points  in  U  U2,  and  for  this  it  suffices  to  consider  the  representation  (1-17) 
for  $*.  Actually,  it  is  more  advantageous  to  work  with  another  representation  which 
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follows  from  (1.17).  For  any  symmetric  tensor  A,  let  A**  denote  the  deviatoric  part 
of  the  square  of  the  deviatoric  part  of  A: 

A"  :=  ((A-)^)-  =  (A-)’  -  14.1  =  (A*)*  -  |4A‘  ,  (1.19) 

where  the  expression  on  the  right  follows  by  squaring  (1.15)  and  taking  the  deviatoric 
part.  Then 

r(A)  =  ^(A)A*+7(A)A** 

=  -i^^.7(A)I  +  ^(A)A*  +7(A)(A*)2,  (1.20) 

where 

^(A)  =  /?(  A)  +  147(A)  .  (1.21) 

The  representation  (1.20)i  is  particularly  useful  in  the  study  of  deviatoric  stress  in 
isotropic  elastic  solids.®  From  (1.21)  it  follows  that  if  0  and  7  are  smooth  on  an  open 
subset  U  of  U  then  so  is  9-,  in  particular,  on  the  isotropic  function  9  is  uniquely 
determined  by  $  and  inherits  any  smoothness  properties  of  However,  it  turns  out 
that  9  is  generally  smoother  than  either  /?  or  7  at  points  in  Ui  U  ZY2.  Indeed,  the  other 
main  result  of  this  paper  is 

Theorem  1.3  Tf  ^  there  is  a  unique  continuous  scalar-valued 

isotropic  function  9  on  U  for  which  the  representations  (1.20)  hold.  If  #  is  (resp. 
C^)  then  9  is  (resp.  C^). 

The  continuity  part  of  this  theorem  is  proved  in  Section  2;  the  differentiability  part 
is  proved  in  Section  4.  From  (1.21)  it  follows  that  /?  is  smooth  if  9  and  7  ^re  smooth; 
we  use  this  fact  and  Theorem  1.3  to  establish  the  smoothness  of  /?  in  Theorems  1.1 
and  1.2. 

The  proofs  of  Theorem  1.2  and  the  second  part  of  Theorem  1.3  use  Proposition  4.1, 
which  gives  conditions  under  which  a  O’"  scalar- valued  isotropic  function  on  can 
be  extended  to  a  function  on  U.  The  proof  of  this  proposition  is  given  in  the 
Section  6;  paxt  of  the  proof  uses  some  results  due  to  Ball  [13].  Section  3  contains 
additional  formulas  relating  9,  7,  and  to  the  first  and  second  derivatives  of  or 
Some  of  these  results  are  used  in  Section  4.  In  the  statements  and  proofs  of 
our  theorems  on  the  smoothness  of  the  coefficients  <1,^,7,  and  9,  we  regard  these 

®For  example,  if  V  is  the  left  stretch  tensor  and  A  =  In  V  is  the  logarithmic  strain  tensor,  then  A* 
and  A**  are  independent  of  the  dilatational  part  J  =  det  V  of  the  deformation.  If  T  =  #(A)  is  the 
Cauchy  stress  tensor,  then  (1.20)i  expresses  the  deviatoric  stress  tensor  as  a  linear  combination  of  the 
volume-independent  deviatoric  strain  tensors  A*  and  A**.  The  coefficients  6  and  7  ma:^be  expr^ed 
as  functions  of  J  (or  lx  —  In  )  and  the  volume-independent,  invariant  strain  measures  and  . 
For  a  hyperelastic  material  with  strain  energy  e  per  unit  reference  volume,  JyiA.)  =  Zde/dlE^ 
and  je{A.)  =  ide/dU^  ,  with  ^(0)  equal  to  twice  the  shear  modulus  of  the  infinitesimal  theory;  cf. 
Scheidler  [11].  An  equivalent  result  in  terms  of  the  representation  (1.20)2  was  obtained  by  Richter 

[12,  (4.6)2]. 
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coefficients  as  (isotropic)  functions  from  U  into  R.  They  may  also  be  regarded  as 
functions  of  the  principal  invariants,  the  moments,  or  the  eigenvalues.  In  Section  5 
we  address  the  question  of  whether  the  established  smoothness  of  the  coefficients, 
regarded  as  isotropic  scalar- valued  functions  on  ZY,  is  inherited  by  the  corresponding 
functions  in  these  representations.  The  answer  depends  on  the  type  of  smoothness 
considered  and  on  which  representation  is  used.  Here  we  also  use  some  results  due  to 
Ball  [13]. 

2  Continuity  of  the  Coefficients 

The  main  results  of  this  section  are  Theorems  2.1-2.3,  which  give  sufficient  conditions 
for  continuity  of  the  coefficients  a,/?, 7,  and  B.  However,  we  begin  with  a  discussion 
of  Serrin’s  example  mentioned  in  the  Introduction.  Define  a  scalar- valued  isotropic 
function  7  on  Sym  by 


7(A)  =  0  if  A*  =0,  7(A)  =  sin(||A*||“^/2)  =  sin(7J]^* if  A*  ^  0  .  (2.1) 

Then  7  is  discontinuous  on  the  set  St/mj  of  all  spherical  tensors,  and  this  discontinuity 
is  nonremovable  since  7  takes  on  all  values  between  —1  and  1  in  every  neighborhood 
of  Sym-^.  Define  a  tensor-valued  isotropic  function  $  on  Sym  by  $(A)  =  7(A)(A*)2. 
Then  $  is  C°°  on  Sym2  U  Sym^,  $  is  differentiable  with  derivative  0  at  each  spherical 
tensor,  and  D^(A)  — »  0  as  A*  —^0;  thus  $  is  on  Sym.  Prom  (1.19)-(1.21)  we 
see  that  ^  =  0  and  /3(A)  =  — |.Za.7(A).  Then  from  (1.18)  and  the  identity 

Va  -  4-  =  jfj ,  (2-2) 

we  have  o:(A)  =  |.Z^7(A).  From  the  uniqueness  of  o:,  ^3,7  on  Sym^  it  follows  that 
these  coefficients  cannot  be  extended  continuously  to  Sym  even  though  $  is  To 
within  some  inessential  constants,  this  is  the  example  given  by  Serrin  [2].  An  even 
more  pathological  example  is  obtained  by  modifying  his  example  as  follows: 

7(A)  =  sin(I^;'^®)  if  A*  ^  0  •  (2.3) 

Then  $  is  still  on  Sym,  but  the  coefficients  a,/3,7  take  on  all  real  values  in  every 
neighborhood  of  Sym^. 

In  the  Introduction  we  observed  that  on  the  open  subset  ZY3  of  ZY,  the  coefficients 
a,  13, 7,  and  B  are  uniquely  determined  by  $  and  inherit  any  smoothness  properties 
of  $.  But  from  the  representations  (1.17)  and  (1.20)  it  is  reasonable  to  expect  that 
on  ZY3  the  coefficients  7,  and  B  are  uniquely  determined  by  and  inherit  any 
smoothness  properties  of^*,  regardless  of  whether  or  not  tr  $  is  smooth.  To  see  that 
this  is  indeed  the  case,  we  can  proceed  as  in  the  derivation  of  (1.14)  and  either  solve 
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(1.17)  for  ^  and  7,  or  solve  (1.20)  for  B  and  7,  and  then  use  the  relation  (1.21).  The 
latter  approach,  together  with  the  identities 

tr  (A* A”)  =  .^* ,  tr  (A**)^  =  ,  (2.4) 

yields 

»(A)  =  [|ff^?tr(«*(A)A-)-3mi.tr($-(A)A")]/Ai.,  (2.5) 

7(A)  =  3[ffi.tr(*’(A)A")-I^.tr(*-(A)A')]/A^.,  (2.6) 

for  A  €  Zfa.  Here  A.,  is  the  cubic  discriminant  of  the  characteristic  polynomial  of 

A*  or,  equivalently,  of  A: 

^A*  =  2^*  -  3^*  =  -  27///]* 

=  K  -  «2)^(a;  -  a;)2K  -  a^)^  =  Aa  ,  (2.7) 

where  a*  axe  the  eigenvalues  of  A*.^° 

If  is  differentiable,  then  9  and  7  are  differentiable  on  74,  so  by  differentiating 
the  representation  (1.20)i  we  obtain  the  following  formula  for  the  derivative  of  at 
any  A  e  74- 

T>$*(A)[E]  =  (V^(A).E)A*  +  (V7(A)*E)A** 

+  ^(A)  E*  +  7(A)(A*E*  +  E*A*)* .  (2.8) 

Here  the  symmetric  tensors  V^(A)  and  V7(A)  are  the  gradients  of  9  and  7  at  A 
(e.g.,  V^(A)  •  E  =  /)^(A)[E]),  and  we  have  used  the  relations 

/)aA*[E]  =  E*,  /\A’*[E]  =  (A*E*+E*A*)*.  (2.9) 

Here  and  below,  E  and  F  denote  arbitrary  symmetric  tensors.  If  is  twice  differ¬ 
entiable,  then  9  and  7  are  twice  differentiable  on  74,  so  by  differentiating  (2.8)  and 
using  (2.9),  we  obtain  the  following  formula  for  the  second  derivative  of  on  74: 

I>2$*(a)[E,F]  =  T>2^(A)[E,F]A*-h/?27(A)[E,F]A** 

+  (V^(A) .  E)  F*  4-  (V7(A)  •  E)(A’F*  4-  F* A*)* 

4-  (V^(A) .  F)  E*  4-  (V7(A)  .  F)(A*E*  4-  E*A*)* 
4-7(A)(F*E*4-E*F*)*.  (2.10) 

We  wish  to  solve  (2.8)  for  ^(A)  and  (2.10)  for  7(A).  The  next  paragraph  contains 
some  algebraic  results  which  will  be  used  in  these  solutions. 

^°Equivalent  formulas  for  6  and  7  were  derived  by  Blinowski  [14]  using  a  slightly  different  proce¬ 
dure;  he  also  discusses  several  representations  which  follow  from  (1.20)i.  We  could  conceivably  use 
(2.5)  and  (2.6)  to  determine  sufficient  conditions  for  the  existence  of  the  limits  of  0(A)  and  7(A)  as 
A  — ^  Wi  U  2/2;  however,  as  indicated  in  the  Introduction,  we  will  follow  a  different  procedure  here. 
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For  any  symmetric  tensor  A,  let  V{A)  denote  the  subspace  of  Sym  consisting  of 

all  polynomials  in  A,  that  is,  all  tensors  of  the  form  CqI  +  c^A  +  C2A?  H - h 

for  some  scalaxs  Cq,  . . . ,  cjv^ .  Then  by  the  Cayley-Hamilton  theorem, 

V{A)  =  span  {I,  A,  A^}  =  span{I,  A*,  A**}  ,  (2-11) 

where  the  right  equality  follows  from  (1.15)  and  (1.19).  Let  'P(A)-*-  denote  the  or¬ 
thogonal  complement  of  ^(A): 

TiA)^  =  {He5ym:I-H  =  A-H  =  A2.H  =  0} 

=  {He5ym:I-H  =  A*-H  =  A**.H  =  0}.  (2.12) 

In  particular,  since  I*H  =  trH,  every  H  €  V{A)-^  is  deviatoric,  i.e.,  H*  =  H.  The 
dimension  of  V{A)  is  ^A,  so  'P(A)-^  has  dimension  6  —  ^A.  In  particular, 

V{A)  =  Sym^  iff  V{Ay  =  Sym*  iff  #A  =  1 .  (2.13) 

For  any  orthonormal  basis  {e,}  =  {61,62,63}  for  V,  let  5'{e,}  denote  the  three- 
dimensional  subspace  of  Sym  consisiting  of  all  symmetric  tensors  with  principal  ba¬ 
sis  {e,};  3 

A  €  5{e,  }  iff  A  =  a,e,-  0  e,- .  (2.14) 

«=i 

For  any  symmetric  tensor  E  let  Eij  =  e,-  •  Ee^-  =  Eji  denote  the  components  of  E 
relative  to  {e^},  where  “  •  ”  also  denotes  the  inner  product  on  V.  Since  a  tensor  H  is 
orthogonal  to  every  A  €  5{eJ  iff  0  =  H  ♦  (e^  0  e,)  =  Hu  for  each  i,  the  orthogonal 
complement  of  5{e,-}  is  the  three-dimensional  subspace  of  symmetric  tensors  whose 
diagonal  components  relative  to  {e,}  axe  zero: 

5{e,}^  =  {H  €  Sym  :  e,-  •  He^  =  Hu  =  0 ,  i  =  1, 2, 3}  .  (2.15) 

In  particular,  each  HG5{ei}-^  is  deviatoric.  If  A  G  S'fe,}  then 

7?(A)c5{e,}  and  5{eJ^  C  7>(A)-l  c  5ym* ,  (2.16) 

with 

7’(A)  =  5{ei}  iff  5{e.}-L  =  7>(A)-L  iff  #A  =  3.  (2.17) 

Also,  for  any  symmetric  tensor  A  we  have 

7^(A)-‘-  =  spanl^le,}-*-  :  {e,}  is  a  principal  basis  for  A} 

=  span{H  :  e,  •  He,  =  0  for  some  principal  basis  {e,}  of  A}.  (2.18) 

The  following  identities  will  be  particularly  useful: 

E*-H  =  E-H  =  tr(EH),  VH  G  5j/m* 

=  VH€S{eiP,  (2.19) 

»<i 
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and 


(A*E*+  E*A*)*. H  =  (A*E  +  EA*)  •  H  =  2tr(A*EH) ,  VH  e  ^(A)-^ 

=  +  a;E,^H^,) ,  V  A  €  5{e.-},  VH  €  5{e.}^ . 

(2.20) 


In  deriving  (2.20)3  have  used  the  fact  that  the  eigenvalues  a*  of  A*  satisfy 


trA*  =  a* +a; +a;  =  0.  (2.21) 

If  is  differentiable,  then  from  (2.8),  (2.12),  (2.19),  and  (2.20),  we  see  that  for 
any  A  €  Us, 

D$*(A)[E]  .H  =  tr(EH)^(A)  +  2tr(A*EH)7(A) ,  VH  G  ViA)^  .  (2.22) 

If  ^  is  an  isotropic  scalar- valued  function  on  U  and  if  U  denotes  the  set  of  all  points  in 
U  at  which  ip  is  differentiable,  then  ZY  is  invariant  and  V'V’  is  an  isotropic  tensor- valued 
function  on  U.  Now  from  the  representation  (1.3)  we  know  that  ^(A)  G  P(A)  for 
any  isotropic  tensor- valued  function  so  by  taking  ^  we  obtain 

V^(A)-K  =  0,  VKg'P(A)-l.  (2.23) 

If  is  twice  differentiable,  then  from  (2.10)  with  F  =  K,  (2.12),  (2.19),  (2.20),  and 
(2.23),  we  see  that  for  any  A  G  U^, 

D2r(A)[E,K].H  =  tr(KH)(V^(A)-E)-f-2tr(A*KH)(V7(A)-E) 

-|-2tr(E*KH)7(A),  VH,K  G  7>(A)-l .  (2.24) 

And  from  (2.24)  and  (2.23),  for  any  A  G  ZY3  we  have 

Z?2$*(A)[L,K]-H  =  2tr(LKH)7(A),  VH,K,L  G  P(A)^  .  (2.25) 

Let  Wle,}  denote  the  subset  of  S'leJ-*-  consisting  of  those  tensors  whose  off-diagonal 
components  relative  to  {e,}  have  absolute  value  1: 


W{e.}  :=  {H  G  Sym  :  Hu  =  0  and  fT.-,-  =  ±1  if  i  /  j}  C  .  (2.26) 

Then  from  (2.19)-(2.21)  with  E  =  H  G  7i{e,},  we  have 

trH2  =  6,  tr(A*H2)  =  0,  V A  G  5{eJ  ,  VH  G  ?f{e.}  .  (2.27) 

On  setting  E  =  H  G  {e,-}  in  (2.22)  and  using  (2.27),  we  obtain  the  following  simple 
formula  for  B  on  Us’- 


B{A)  =  |Pr(A)[H] .  H ,  VH  G  W{e.} ,  (2.28) 
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where  {ej  is  any  principal  basis  for  A.  Similarly, 

trH3  =  6£H,  £«  :=  ,  VHe5{e,}-L,  (2.29) 

so  on  setting  L  =  K  =  H  G  H{e,}  in  (2.25)  and  using  (2.29)  and  (2.26),  we  obtain 
the  following  simple  formula  for  7  onU^: 

7(A)  =  gfl2«-(A)lH,H]-H,  £h  =  ±1.  VH  €«{£;},  (2.30) 

where  {e,}  is  any  principal  basis  for  A. 

K  is  then  is  continuous,  and  we  can  use  (2.28)  to  show  that  ^(A)  has 
a  limit  as  A  — »  B  G  U  2/2  from  within  U^.  The  only  difficulty  is  that  as  A  — >  B 
along  some  arbitrary  path  in  ZY3,  the  principal  axes  of  A  will  generally  change;  in 
particular,  they  might  not  have  limits  at  B.^^  Since  H  in  the  formula  (2.28)  belongs 
to  the  set  TfleJ,  which  changes  with  the  principal  basis  {e,}  of  A,  it  is  not  obvious 
that  (2.28)  has  a  limit  at  each  point  B  G  U  ZY2.  Similar  comments  apply  to  the 
formula  (2.30)  for  7  when  is  Fortunately,  this  difficulty  is  easily  overcome  by 
utilizing  the  isotropy  of  B  and  7,  as  we  now  show. 

For  any  orthonormal  basis  {e,}  for  V,  the  set  fl  5{e,}  consists  of  all  tensors 
in  U  with  three  distinct  eigenvalues  and  principal  basis  {e^}.  Included  in  this  set  is 
the  set  ZY<{e,}  consisting  of  all  A  =  a,e,-  (g)  e,-  G  ZY  for  which  Cj  <  02  <  G3.  The 

invariance  of  ^^3  guarantees  that  ZY3  fl  5{e,}  and  are  nonempty.  And  since  ZY3 

is  an  open  subset  of  the  six-dimensional  space  Sym,  lAz  D  5{e,}  and  ZY^{e,}  are  open 
subsets  of  the  three-dimensional  subspace  5{e,}. 

Proposition  2.1  For  a  continuous  scalar-valued  isotropic  function  ip  onli^,  the 
following  conditions  are  equivalent: 

(1)  Ip  has  a  continuous  extension  ip  to  U; 

(2)  VB  G  ZYi  U  U2,  the  lim  ^(A)  exists; 

A->B 

AeWs 

(3)  VB  G  U  Z/2  and  each  principal  basis  {e,}  for  B,  the  lim  ip{A)  exists; 

A-i-B 

A€W3n5{e,} 

(4)  VB  G  ZYi  U  Z^2  and  some  principal  basis  {61,62,63}  for  B  corresponding  to  the 

ordered  eigenvalues  b^  <  b2  <  b^  of  B,  the  lim  ip{A)  exists. 

A->B 

A€W<{e.} 

When  these  conditions  hold,  ip  is  unique  and  isotropic  on  U. 

^^Cf.  the  example  in  Kato  [15,  p.  128]  or  the  two  examples  in  Scheidler  [16]. 
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As  discussed  in  the  Introduction,  the  equivalence  of  (1)  and  (2)  follows  from  the 
fact  that  is  dense  in  i/,  and  when  these  conditions  hold  ^  is  unique:  ^(B)  = 
limA._^B  V’(-A-),  A  e  Z4,  for  each  B  G  U  U2.  Clearly,  (2)  ^  (3)  =4>  (4),  so  it  remains 
to  show  that  (4)  implies  (2)  and  that  the  continuous  extension  ^  is  isotropic.  The 
proof  will  be  completed  below, f[j.st  we  apply  Proposition  2.1  to  the  continuity 
of  B  and  7.  Condition  (3)  will  suffice  for  all  the  applications  in  this  paper.  However, 
the  weaker  condition  (4)  arises  naturally  in  the  proof  below  and  has  other  applications 
that  axe  not  discussed  here. 

Let  B  G  Wj  U  ZY2  and  let  {e,}  be  any  principal  basis  for  B.  Choose  any  H  G 
7f{ei};  then  the  formula  for  ^(A)  in  (2.28)  holds  for  each  A  G  W3  D  <S'{e,}.  If  #* 
is  then  is  continuous,  so  the  right-hand  side  of  (2.28)  has  a  limit,  namely 
|D#*(B)[H]  •  H,  as  A  ^  B  from  within  fl  ^{e,},  since  H  is  fixed.  Thus  condition 
(3)  of  Proposition  2.1  holds  for  ^  and  we  know  that  B  is  continuous  and  isotropic 
on  Uzi  so  by  Proposition  2.1  we  conclude  that  B  has  a  unique  continuous  extension  to 
U.  Summarizing  (and  replacing  B  by  A  in  the  above  limit),  we  have 

Theorem  2.1  If  is  on  U,  then  there  is  a  unique  continuous  scalar-valued 
isotropic  function  B  on  U  for  which  the  representations  (1.20)  hold.  This  coefficient 
B  is  given  by  (2.28)  for  any  A  Gld  and  any  principal  basis  {e,}  for  A. 

By  a  similax  axgument  applied  to  (2.30),  we  conclude  that  7  has  a  continuous  exten¬ 
sion  to  U  when  is  C^.  On  combining  this  result  with  Theorem  2.1  and  using  the 
relations  (1.18)  and  (1.21),  we  obtain 

Theorem  2.2  If  is  on  Li,  then  there  are  unique  continuous  scalar-valued 
isotropic  functions  l3,‘y,B  onli  for  which  the  representations  (1.3),  (1.17),  and  (1.20) 
hold.  The  coefficients  7  and  B  are  given  by  (2.30)  and  (2.28),  respectively,  for  any 
A  £11  and  any  principal  basis  {e,}  for  A;  then  ^  can  be  expressed  in  terms  of  B  and 
7  by  (1.21).  If  in  addition,  tr  #  is  continuous,  then  there  is  a  unique  continuous 
scalar-valued  isotropic  function  a  onli  for  which  the  representation  (1.3)  holds,  and 
a  is  given  in  terms  of  tr$,  and  7  by  (1.18). 

Proof  of  Proposition  2.1:  Assume  (4)  holds,  so  that  B  =  6,e,-  (g>  e,-  with 

b\  <b2<  63.  Let  /b  denote  the  hmit  in  (4).  Let  {A„}  (n  =  1,2,.. .)  be  any  sequence 
in  Z/3  converging  to  B,  and  let  ai,„e,-^n0e,-^„,  with  <  02, „  <  03, „,  be  a  spectral 
decomposition  of  A„;  then  a,-  „  6^.^^  Let  Q„  be  the  orthogonal  tensor  which  maps 

e,-^„  to  e,-  {i  -  1,2,3).  Then  Q„A„Q^  =  a.-,ne,  ®e.-  G  24<{e,},  and  since  6.- 

it  follows  that  Q„A„Q^  — >  B.  Thus  (4)  implies  that  V’(QnA„Q^)  ^  /b-  But  then 

^^Proposition  2.1  also  follows  from  Proposition  6.1.  We  have  given  a  separate  proof  of  Proposi¬ 
tion  2.1  in  this  section  to  avoid  the  discussion  of  certain  preliminary  results  which  are  only  needed 
for  the  more  general  Proposition  6.1. 

^^Here  we  have  used  the  well-known  fact  that  the  ordered  eigenvalues  of  a  symmetric  tensor  are 
continuous  functions  on  Sym,  In  fact,  they  are  Lipschitz  continuous;  a  proof  can  be  found  in  Ball 
[13,  Lemma  5.8]. 
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V’(A„)  — ^  /b  since  ^(A„)  =  V’(QnA„Q]|)  by  the  isotropy  of  Since  {A„}  was  an 
arbitrary  sequence  in  converging  to  B,  the  limA_,B  V’(A),  A  G  Z4,  exists  and  equals 
Zb-  In  particular,  (2)  holds. 

Since  the  continuous  extension  •0  coincides  with  the  isotropic  function  ^  on  ZY3, 
to  prove  the  isotropy  of  ^  on  ZZ  it  suffices  to  show  that  is  isotropic  on  ZYi  U  ZY2-  Let 
B  G  ZYi  U  U2-,  let  Q  be  any  orthogonal  tensor,  and  let  {B„}  be  any  sequence  in  Z4 
converging  to  QBQ^.  Then  'i/’(B„)  — >  ^(QBQ^),  but  also  ^(B„)  =  ^(Q^B„Q)  — > 
^(B),  since  Q'^B^Q  B.  Thus  ^^(QBQ^)  =  □ 

If  $  is  C*"  then  so  are  and  tr#.  Therefore,  from  Theorem  2.1  we  recover  the 
first  statement  in  Theorem  1.3,  and  from  Theorem  2.2  we  recover  the  case  where  $ 
is  on  H  in  Theorem  l.!.^"*  The  next  theorem  implies  the  case  where  #  is  on 
U2  U  Z4  in  Theorem  1.1.  Note  that  we  cannot  use  the  formula  (2.30)  for  7  in  this 
case  since  we  axe  not  assuming  that  is  twice  differentiable.  Let  B  =  Z»,e,-  0  e,- 

be  any  spectral  decomposition  of  B  G  ZY2.  Then  bi  ^  bj„  =  b„  for  some  permutation 
l,m,n  of  1,2,3.  Let  H,K  G  satisfy 


Ki^  =  0,  =  =  =  =  (2.31) 

with  Him  arbitrary.  Then  from  (2.19)  and  (2.20),  we  have 

tr(KH)  =  0  and  tr(A*KH)  =  a;  -  =  a,  -  c,„  ,  V  A  G  5{e,}  .  (2.32) 

On  setting  E  =  K  in  (2.22),  we  find  that  for  any  A  G  Z^3  n  5{e,}, 

T>$*(A)[K]-H 

2(a/  <^tn) 

Now  as  A  ^  B  from  within  Z4  H  5'{e,},  O;  =  e/- Ae^  —>■  e/*Be/  =  ;  similarly, 

Om  Since  bi  —  b^  ^  0,  and  since  H  and  K  may  be  held  fixed  during  this 

limiting  process,  we  see  that  limA_B7(A),  A  G  Z4  D  5'{e,},  exists  when  Z)$*  is 
continuous.  And  we  know  that  7  is  continuous  and  isotropic  on  ZY3.  By  applying 
Proposition  2.1  with  ^  =  7  and  ZZ  replaced  with  U  :=U  —  Ui  =  ZY2  U  Z^  (so  that 
Ux  =  0,  Z4  =  ^2?  ^3  =  ^3)?  we  see  that  7  has  a  continuous  extension  to  ZY2  U  ZZ3.  On 
combining  this  result  with  Theorem  2.1  and  replacing  B  by  A  in  the  above  limit,  we 
obtain 

^^The  weaker  conditions  on  #  in  Theorems  2.1  and  2.2  (i.e.,  the  requirement  that  #*  is 
while  tr  $  need  only  be  continuous)  may  have  applications  to  isotropic  elastic  solids.  Suppose,  for 
example,  that  the  material  has  suffered  damage  in  the  form  of  a  random  distribution  of  microcreicks. 
Assuming  that  the  damaged  material  can  also  be  modeled  as  an  isotropic  elastic  solid,  one  expects 
a  lower  bulk  modulus  in  tension  (where  cracks  are  opening),  while  in  compression  (where  all  cracks 
are  closed)  the  bulk  modulus  should  ideally  be  the  same  as  in  the  undamaged  material;  cf.  Horii  and 
Nemat-Nasser  [17],  where  an  analysis  based  on  the  linear  theory  also  predicts  that  the  shear  modulus 
in  the  damaged  material  increases  smoothly  in  going  from  tension  to  compression.  Thus  if  $  is  the 
response  function  for  the  Cauchy  stress  tensor  in  the  damaged  material,  and  hence  —  |tr#  and 
the  response  functions  for  the  pressure  and  deviatoric  stress,  respectively,  it  might  be  reasonable  to 
regard  tr  $  as  continuous  but  with  a  discontinuous  derivative,  and  #*  as  (fc  >  1). 


,  VH,K  G  5{e.}-L  satisfying  (2.31).  (2.33) 
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Theorem  2.3  If  is  on  the  open  set  U  —  Ux  —  U2U  U^,  then  there  are 
unique  continuous  scalar-valued  isotropic  functions  /3, 7, 6  on  112  U  ^3  which  the 
representations  (1.3),  (1.17),  and  (1.20)  hold.  The  coefficients  7  and  9  are  given  by 
(2.33)  and  (2.28),  respectively,  for  any  A  G  i/2  U  Z/3  and  any  principal  basis  {e,}  for 
A;  in  particular,  the  ai  and  a^  in  (2.33)  are  the  distinct  eigenvalues  of  A.  Then  ^ 
can  be  expressed  in  terms  of  9  and  7  by  (1.21).  If,  in  addition,  tr$  is  continuous 
on  U2  U  U2,  then  there  is  a  unique  continuous  scalar-valued  isotropic  function  a  on 
U2  U  i/3  for  which  the  representation  (1.3)  holds,  and  a  is  given  in  terms  of  tr#, 
and  7  by  (1.18). 

3  Alternate  Formulas 

Here  we  consider  alternate  formulas  for  9  and  7  and  for  the  derivatives  of  $  and  $*. 
Some  of  these  results  are  used  in  the  next  section.  Throughout  this  section,  E,  F,  G 
denote  arbitrary  symmetric  tensors,  and  l,m,n  denotes  an  arbitrary  permutation  of 
1,2,3.  For  any  orthonormal  basis  {e,}  =  {©1,62,63}  for  V,  let 

H,  :=  e„,  0  e„  +  e„  0  e„,  e  5{e,}-L .  (3.1) 

Then 

H,-H,=2,  H,-H„  =  0;  (3.2) 

hence  {Hi,H2,H3}  is  an  orthogonal  basis  for  S'fe, Also, 

=  tr(HiH„.H„)  =  l,  tr(HiH2)  =  0.  (3.3) 

Since  the  eigenvalues  of  A*  satisfy  (2.21),  for  any  A  €  5{e,}  we  have 

A*  Hi  +  H,A*  =  -a;Hj ,  (3.4) 

(A*Hi  +  H,A*)  •  Hi  =  2  tr(A*H2)  =  -2a;  ,  (3.5) 

and 

(A*Hi  +  HiA*)  •  H,„  =  2  tr(A*HiH„.)  =  0 .  (3.6) 

Referring  to  the  definition  (2.26)  of  7/{e,  },  we  see  that 

Hi  +  H2  +  H36W{eJ.  (3.7) 

Assume  is  differentiable,  and  let  A  €  i/3.  From  (2.8)  and  (2.23)  we  obtain 

T>$*(A)[H]  =  ^(A)H  +  7(A)(A*H4-HA*),  VH  g  P(A)-^,  (3.8) 

where  we  have  used  the  fact  that  H  and  A*H  +  HA*  are  deviatoric.  We  can  use 
(1.15)  and  (1.21)  to  express  I?$*(A)[H]  in  terms  of  /?  instead  of  9: 

Z)#*(A)[H]  =  ^(A)H  +  7(A)(AH  +  HA),  VH  g  ViA)^.  (3.9) 
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Similaxly,  from  (2.22),  (1.15),  and  (1.21),  we  have 

I>r(A)[E]  •  H  =  tr(EH)/?(A)  +  2  tr(AEH)7(A) ,  V H  e  ^(A)-^.  (3.10) 

On  taking  the  inner  product  of  (3.8)  or  (3.9)  with  E,  we  obtain  the  same  expression 
as  on  the  right-hand  side  of  (2.22)  or  (3.10),  respectively.  Therefore 

I>$*(A)[H]-E  =  !)$*(  A)[E]-H,  VH6  7?(A)-^.  (3.11) 

Since  A*C  and  AC  €  V{A)  for  any  C  €  7?(A),  by  setting  E  =  C  in  (2.22)  or  (3.10) 
and  using  (3.11),  we  obtain 

T>$*(A)[H].C  =  I>$*(A)[C]-H  =  0,  VC  G  P(A),  VH  G  7?(A)-L.  (3.12) 

By  setting  E  =  H/  and  H  =  H„,  in  (2.22)  and  using  (3.2)  and  (3.6),  we  obtain 

Zir(A)[Hi].H„,:=0.  (3.13) 

By  (3.7)  we  may  set  H  =  in  (2.28),  and  when  (3.13)  is  applied  to  this  result, 

we  obtain 

«(A)  =  i  E  CS'CAjlHJ  ■  Ht .  (3.14) 

k=\ 

In  (3.14)  and  (3.12),  and  in  similar  results  below,  the  orthonormal  basis  {e^}  relative 
to  which  Hi,H2,  H3  are  defined  is  assumed  to  be  a  principal  basis  for  A.  Also,  since 
only  differentiablity  of  was  assumed  above,  the  relations  (3.8)-(3.14)  are  generally 
valid  only  for  A  G  7/3. 

Now  assume  that  $  is  differentiable.  From  (1.16)  we  have 

P$(A)[E]  =  I>$*(A)[E]  +  |[V(tr#)(A)  -Ejl.  (3.15) 

Since  I  •  G*  =  0  for  any  tensor  G,  it  follows  that 

I>$(A)[E] .  G*  =  £»#*(A)[E]  •  G*  .  (3.16) 

From  (3.15)  and  (2.23)  with  ^  =  tr#,  it  follows  that 

I>#(A)[H]  =  !>$*( A)[H] ,  V H  G  V{A)^  .  (3.17) 

These  results  hold  for  each  A  eU.  Now  assume  that  A  G  Then  by  (3.17)  we  may 
replace  by  $  in  (3.8)  and  (3.9).  Similarly,  since  H*  =  H  for  each  H  G  V{A)-^, 
by  (3.16)  and  (3.17)  we  may  replace  by  $  in  (2.22),  (2.28),  (2.33),  and  (3.10)- 
(3.14).  K  a  scalar-valued  isotropic  function  V’  on  7/  is  twice-differentiable  on  T/j,  then 
by  applying  (3.12)  and  (3.13)  with  =  VV’,  we  see  that  for  any  A  G  7/3, 

I>V(A)[C,H]  =  0,  VC  G  7>(A),  VH  g  7?(A)^,  (3.18) 
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and 


(3.19) 


C“V’(A)[H„  H„]  =  0 ,  VV-(A)  ■  H,  =  0 , 

where  the  relation  on  the  right  is  just  a  special  case  of  (2.23). 

Now  assume  that  is  C^.  Then  is  continuous,  and  without  loss  of  generality 
we  may  assume  that  6  is  also  continuous  on  U  (cf.  Theorem  2.1).  On  taking  the  limit 
of  (3.13)  or  (3.14)  as  A  — ^  B  =  ®  e,-  G  Z/j  U  ZY2  from  within  Uz  fl  5{e,}  while 

holding  the  fixed,  and  then  writing  A  in  place  of  B,  we  see  that  (3.13)  and  (3.14) 
hold  for  every  A  €  U.  By  applying  the  same  limiting  process  to  (3.11)  and  (3.12) 
and  using  (2.17),  we  see  that  for  every  A  eU  with  principal  basis  {e,}, 

I>r(A)[H]  .C  =  Z1$*(A)[C]  .H  =  0,  VC  e  5{e.},  VH  €  5{e,}J-,  (3.20) 

and  D$*(A)[H]*E  =  i}$*(A)[E]*H  for  each  H  G  5{e,}-^.  In  the  latter  relation 
E  is  arbitrary,  so  from  (2.18)  we  conclude  that  (3.11)  holds  for  each  A  Eli.  From 
the  previous  paragraph  it  follows  that  these  results  also  hold  with  replaced  by  #, 
provided  that  $  is  or  more  generally,  that  is  and  tr$  is  differentiable, 

since  the  latter  conditions  suffice  for  (3.16)  and  (3.17). 

Now  assume  that  is  twice  differentiable,  and  let  A  G  Uz-  By  setting  E  =  H/ 
and  F  =  in  (2.10)  and  using  (3.19)  and  (3.3)i,  we  obtain 

D^^\A)[Hi,  H„]  =  7(A)  H„  .  (3.21) 

On  taking  the  inner  product  of  (3.21)  with  H„  or  and  using  (3.2),  or  from  (2.25) 
and  (3.3),  we  obtain 

7(A)  =  (3.22) 

and 

D2*-(A)[H„H„]-H„=0.  (3.23) 

Now  assume  that  #  is  twice  differentiable.  Then  from  (3.15)  we  have 

P2$(A)[E,  F]  =  I>2$*(a)[E,  F]  +  |Z)2(tr  $)(A)[E,  F])  I ,  (3.24) 

which  implies 

D2$(a)[E,  F]  .  G*  =  D2$*(A)[E,  F] .  G*  .  (3.25) 

These  results  hold  for  any  A  G  7/.  Now  assume  A  G  Hz-  Then  from  (3.25)  it  follows 
that  we  may  replace  by  $  in  the  formulas  (2.24),  (2.25),  (2.30),  (3.22),  and  (3.23). 
By  setting  E  =  H/  and  F  =  in  (3.24)  and  using  (3.19)  with  V’  =  tr$,  we  obtain 

r>2*{A)[H„  H„1  =  i)2r(A)[H,,  H„] .  (3.26) 

Thus  (3.21)  holds  with  replaced  by  $. 

Now  assume  that  is  C^.  Then  is  continuous,  and  without  loss  of  gen¬ 

erality  we  may  assume  that  0,  7,  and  ^  axe  also  continuous  on  U  (cf.  Theorem  2.2). 
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Arguing  as  in  the  case  above,  we  conclude  that  (3.8)-(3.10)  and  (3.21)-(3.23)  hold 
for  every  A  G  And  from  the  previous  paragraph  it  follows  that  these  results  hold 
with  $*  replaced  by  $  when  #  is  C^.  Actually,  (3.8)-(3.10)  with  $  hold  for 

every  A.  Eli  under  the  weaker  condition  that  is  and  tr  $  is  differentiable,  since 
this  condition  suffices  for  (3.16)  and  (3.17).  Similarly,  (3.21)-(3.23)  with  $ 

hold  for  every  A  E  U  under  the  weaker  condition  that  is  and  tr$  is  twice 
differentiable. 


4  Differentiability  of  the  Coefficients 

To  establish  sufficient  conditions  for  the  differentiablity  of  the  coefficients  we  use  the 
same  approach  as  in  the  previous  section,  that  is,  for  A  G  Z/3  we  express  V0(A)  and 
V7(A)  in  terms  of  Z)2#*(A)  or  Z?^#*(A)  and  show  that  these  expressions  have  limits 
as  A  — >  B  G  T/j  U  7/2  from  within  D  S'fe,},  where  {e,}  is  a  principal  basis  for 
B.  Then  to  conclude  that  0  and  7  are  functions  on  ZY,  we  need  an  appropriate 
generalization  of  Proposition  2.1.  Since  we  are  also  interested  in  analogous  results 
for  higher  order  derivatives  of  the  coefficients,  we  generalize  Proposition  2.1  to  C’" 
isotropic  functions,  where  r  denotes  some  positive  integer. 

Proposition  4.1  Let  ip  be  a  scalar-valued  isotropic  function  onlA  which  is 
C’'  on  ZY3.  Then  the  following  conditions  are  equivalent: 

(1)  Ip  is  C’’  on  U; 

72)  VB  G  ZYj  U  ZY2,  the  lim  D’'ip{A)  exists; 

A->B 

Aei4 

(3)  VB  G  ZYi  U  U2,  each  principal  basis  {ej  for  B,  and  any  Ei, . . . , G  Sym, 

the  lim  £)’'^(A)[Ei, . . . ,  E^]  exists. 

A-fB 

Ae%nS{e<} 

Here  (7°  (i.e.,  for  r  =  1)  means  continuous.  The  proof  is  given  in  Section  6.^® 
Throughout  this  section  we  assume  that  is  at  least  on  ZY.  Then  without  loss  of 
generality  we  can  assume  that  7,  and  6  are  continuous  on  ZY  (cf.  Theorem  2.2).  As 
in  the  previous  section,  E,  F,  G  denote  arbitrary  symmetric  tensors.  Unless  stated 
otherwise,  /,  m,  n  denotes  an  arbitrary  permutation  of  1,2,3. 

By  setting  K  =  H  G  7Y{e,}  in  (2.24)  and  using  (2.27),  we  obtain  the  following 
formula  for  V0  on  ZY3: 


V^(A)-E  =  |D2$*(a)[E,H]-H- itr(E*H2)7(A),  VHG?Y{e.},  (4.1) 

^®As  shown  in  Section  6,  condition  (3)  can  be  replaced  by  a  much  weaker  condition.  However,  the 
version  stated  above  suffices  for  the  applications  in  this  paper. 
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where  {e,}  is  any  principal  basis  for  A.  Let  B  €  ZYj  U  i/2  and  let  {e,}  be  any  principal 
basis  for  B.  Choose  any  H  €  W{e,}  and  any  symmetric  E.  Then  the  formula  for 
V^(A)  in  (4.1)  holds  for  every  A  €  i/3  fl  5{e,}.  Since  D^*  and  7  are  continuous  on 
i/,  the  right-hand  side  of  (4.1)  has  a  limit  as  A  — »■  B  from  within  i/3  fl  5{e,},  and 
thus  so  does  D^(A)[E]  =  V^(A)  ^E.  Since  9  is  isotropic  and  continuous  on  U  and 

(in  fact  C^)  on  i/3,  Proposition  4.1  with  xj)  =  9  and  r  =  1  yields' 

Theorem  4.1  If  is  on  U,  then  9  is  on  U,  and  V9  is  given  by  (4-1)  for 
any  A  G  i/  and  any  principal  basis  {e,}  for  A. 

Now  assume  that  is  C^.  By  differentiating  (2.10)  and  using  (2.9),  we  obtain 
the  following  formula  for  D^#*(A)  at  any  A  €  i/3: 

Z?3$*(A)[E,  F,  G]  =  i:>3^(A)[E,  F,  G]  A*  I>37(A)[E,  F,  G]  A** 

-h  L>2^(A)[E,F]  G*  +  £»27(A)[E,F](A*G*  -f  G*A*)* 

-1-  D^9{A)[E,  G]  F*  -f  L>2^(A)[E,  G](A*F*  -1-  F*A*)* 

-f-  I>2^(A)[F,  G]  E*  -f  D^'r{A)[F,  G](A*E*  +  E*A*)* 
d-  (V7(A)  .  E)(G*F*  -t-  F*G*)* 

(V7(A)  •  F)(G*E*  -f  E*G*)* 

-H  (V7(A)  •  G)(E*F*  +  F*E*)* .  (4.2) 

On  taking  the  inner  product  of  (4.2)  with  H  G  V{A)^  and  using  (2.12),  (2.19),  and 
(2.20),  we  obtain 

D3#*(a)[E,F,G].H  =  D2^(A)[E,F]tr(GH)  +  2£>27(A)[E,F]tr(A*GH) 

I>2^(A)[E,  G]  tr(FH)  +  2D^  j{A)[E,  G]  tr(A*FH) 

+  L>2^(A)[F,  G]  tr(EH)  -f  2I>2^(A)[F,  G]  tr(A*EH) 

-h  2(V7(A)  •  G)  tr(E*F*H)  -h  2(V7(A)  •  F)  tr(G’E*H) 

+  2(V7(A)-E)tr(G*F*H),  (4.3) 

for  any  A  G  i/3  and  any  H  G  ^(A)-*-.  We  wish  to  solve  this  for  V7(A)*E.  Set 
F  =  H/,  G  =  Hto,  and  H  =  H„.  Then  the  first  two  lines  in  (4.3)  drop  out  from  (3.2) 
and  (3.6);  and  the  third  and  fourth  lines  drop  out  from  (3.19).  On  using  (3.3)2  in  the 
last  line,  we  obtain  the  following  simple  formula  for  V7  on  i/3: 

V7(A)  .  E  =  iD3#*(A)[E,  Hi,  H,„]  ♦  H„ ,  (4.4) 

where  the  orthonormal  basis  {e,-}  relative  to  which  Hi,H2,H3  are  defined  (cf.  (3.1)) 
is  any  principal  basis  for  A.  Arguing  as  in  the  proof  of  Theorem  4.1,  we  conclude 
that  7  is  on  14  when  is  and  that  (4.4)  holds  for  each  A  G  i/. 

Next  we  wish  to  show  that  9  is  on  i/  when  is  C^.  Since  9  is  on  U 
(cf.  Theorem  4.1),  by  Proposition  4.1  we  need  only  show  that  limA_.B  D‘^9{A)\E,  F], 
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A  €  Z4  n  5{e,},  exists  for  each  B  €  ZYi  U  Z/j,  each  principal  basis  {e,}  for  B,  and  any 
symmetric  tensors  E,F.  And  since  {Hi,H2,H3}  is  a  basis  for  it  suffices  to 

show  that  D2^(A)[H,,HJ,  D^e{A)[Hi,Hi],  Z>2^(A)[C,D],  and  I?2^(A)[C,Hi]  have 
limits  from  within  2/3  D  S'leJ,  for  any  C,D  G  5'{e,}.  But  the  first  and  last  of  these 
expressions  are  zero  by  (3.19)  and  (3.18),  respectively.  Consider  the  third  expression. 
Set  E  =  C,  F  =  D,  and  G  =  K  in  (4.3),  where  C,  D  G  V{A)  and  K  G  V{A)^.  Then 
A*C,A*D,  C*D*  G  'P(A),  so  the  second  and  third  lines  and  the  first  expression  in 
the  fourth  line  in  (4.3)  drop  out  since  H  G  V{A)-^.  Thus  for  any  A  G  W3  we  have 

I>3$*(A)[C,D,K].H  =  Z)2^(A)[C,D]tr(KH)  +  2D27(A)[C,D]tr(A*KH) 

+  2(V7(A)  .  D)  tr(C*KH)  +  2(V7(A)  •  C)  tr(D*KH) , 
VC,Dg7>(A),  VH,Kg7?(A)^.  (4.5) 

Now  let  K  =  H  G  'H{e,}  in  (4.5),  where  {e,}  is  any  principal  basis  for  A  G  U3.  Then 
V{A)  =  S'fe,},  so  from  (2.27)  we  obtain 

£>2^(A)[C,D]  =  iZl3$«(A)[C,D,H].H,  VC,D  G  5{e.},  VH  G -^{e,} .  (4.6) 

Since  is  it  follows  that  has  a  limit  as  A  ^  B  G  U  7^2 

from  within  U3  D  5{e;},  with  C,D  G  5{e,}  fixed.  Finally,  consider  the  expression 
I?20(A)[H/,  H/].  On  settting  E  =  F  =  and  G  =  H  =  H/  or  H„,  in  (4.3),  and 
using  (3.2),  (3.5),  and  (3.19),  we  obtain 

1D3#*(A)[H,,  H,] .  H,  =  D^e{A)[lli,  -  a*Zl2^(A)[H„  H,]  (4.7) 

and 

1Z)3«-(A)[H,.  H„  HJ  •  H„  =  I>=«(A)1H„  H,J  -  o;.D=7(A)[H„  H,] .  (4.8) 

On  axiding  (4.7)  and  (4.8)  and  using  —a^  °n’  obtain 

2D2«(A)(H„H,]+o;£|27(A)[H„H,]  =  4D“«*(A)[H,.H,,H,].H,+ 

lD=r(A)|H,,H,,HJ-H„.  (4.9) 

Then  by  interchanging  m  and  n  in  (4.9)  and  adding  the  result  to  (4.8),  we  obtain 

3D2^(A)[H„H,]  =  lI>5#*(A)[H„H„H„.].H,„  +  p3$*(A)[H„H„H„].H„ 

+  |£»3$*(a)[H„H„H,]-H,,  (4.10) 

where  the  orthonormal  basis  {e,}  relative  to  which  Hi,H2,H3  are  defined  is  any 
principal  basis  for  A.  It  follows  that  £)2^^A)[H;,  H;]  has  a  limit  as  A  — »  B  G  ZYj  U  7^2 
from  within  74  fl  5{e,}  with  H/  fixed.  From  the  remarks  at  the  beginning  of  the 
paragraph,  we  conclude  that  6  is  on  U  when  is  C^.  On  combining  this  with 
the  result  of  the  previous  paragraph  and  using  (1.21)  and  (1.18),  we  obtain 
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Theorem  4.2  If  is  on  U,  then  0  is  on  U,  and  7  and  ^  are  on  Id. 
For  any  A  gU,  Vd(A)  and  V7(A)  are  given  by  (4.I)  and  (4-4))  respectively,  and 

V/?(A) .  E  =  V^( A) .  E  -  |(  Je7( A)  +  /aV7(A)  •  E) .  (4.11) 

If,  in  addition,  tr$  is  on  U,  then  so  is  a,  and  Va  can  be  expressed  in  terms  of 
V(tr$),  VyS,  and  V7  by  means  of  the  relation 

V(tr$)(A)-E  =  (3Va(A)  +  4V^(A)  +  irAV7(A))  •  E 

+  /b^(A)  +  2(A.E)7(A).  (4.12) 

Now  we  turn  to  the  smoothness  of  7  at  points  in  142-  Let  B  =  &,e,'  0e,-  be  any 

spectral  decomposition  of  B  €  142-  Then  bi  ^  b^  =  b„  for  some  permutation  I,  m,  n  of 
1,2,3.  For  H,K  G  *S'{e,}-L  satisfying  (2.31)  we  have,  by  (2.32)  and  (2.24), 

2{ai  -  aJV^iA)  •  E  =  i}2$*(A)[E,K]  -  H  -  2tr(E*KH)7(A)  (4.13) 

for  any  A  £  14^  with  principal  basis  {e,}.  Then  from  (4.13),  Proposition  4.1,  and 
arguments  similar  to  those  preceding  Theorem  2.3,  we  conclude  that  7  is  on 
142  L*  14^  if  is  on  142  U  2^3-  Now  assume  that  is  on  i/2  U  i/3.  Then  V7  is 
continuous  on  i/2  U  i/3,  and  for  H  and  K  as  above,  (4.5)  yields 

T»3r(A)[C,D,K]-H  =  2(a,-a„,)D27(A)[C,D]  +  2(q-c,„)(V7(A).D) 

+  2(d,-4)(V7(A).C)  (4.14) 

for  any  A  e  i/3  D  5{eJ  and  any  C,D  €  ^{ej.  By  subtracting  (4.7)  from  (4.8),  we 
obtain 

(a;-a„,)£>27(A)[Hi,H,]  =  iT>3#*(A)[H„H„Hj  •  H,„ 

-iI?3r(A)[Hi,H„H,].H,.  (4.15) 

Arguing  as  above  (see  also  the  comments  preceding  (4.5)),  we  conclude  that  7  is 
on  i/2  U  I4z-  On  combining  the  results  of  this  paragraph  with  Theorems  4.1  and  4.2 
and  using  (1.21)  and  (1.18),  we  obtain 

Theorem  4.3  If  is  (resp.  C^)  on  i/2  U  i/3,  then  ^,7,jS  are  (resp.  C^) 
on  i/2  U  i/3.  If,  in  addition,  tr#  is  (resp.  C'^)  on  U2  U  14^,  then  so  is  a. 

Of  course,  if  $  is  C*  then  and  tr$  are  C^,  so  Theorems  4.1,  4.2,  and  4.3 
yield  Theorem  1.2  and  the  second  part  of  Theorem  1.3.  From  (3.25)  it  follows  that 
we  may  replace  with  $  in  (4.1)  and  (4.13)  when  $  is  C^.  If  is  C^,  then  by 
differentiating  (3.24)  we  obtain 

Z)3$(A)[E,  F,  G] .  L*  =  T»3$*(A)[E,  F,  G]  •  L*  (4.16) 

for  any  symmetric  tensors  E,  F,  G,L,  so  we  may  replace  by  #  in  (4.3)-(4.10), 
(4.14),  and  (4.15). 
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5  Smoothness  of  the  Representations  for  the  Coef¬ 
ficients 

If  V’  is  a  scalar- valued  isotropic  function  on  U,  then  there  are  scalar- valued  functions 
•if),  if),  if),  and  on  appropriate  subsets  of  R^,  such  that 

V.(A)  =  = 

=  <A(oi,<i2.<i3)  =  i(’.(4.J4*>-^-f4').  (5-1) 

with  if)  symmetric  in  its  arguments.^®  Conversely,  li  •i})  :U  U  has  any  one  of  the 
representations  above,  then  if)  is  isotropic.  In  the  statements  and  proofs  of  our  the¬ 
orems  on  the  smoothness  of  the  coefficients  0,^,7,  and  6,  we  have  regarded  these 
coefficients  as  functions  on  U.  Of  course,  since  the  coefficient  functions  are  isotropic, 
they  each  have  representations  of  the  form  (5.1).  We  are  then  faced  with  the  ques¬ 
tion  of  whether  the  established  smoothness  of  the  coefficients,  regarded  as  isotropic 
scalar- valued  functions  on  U,  is  inherited  by  the  corresponding  functions  in  the  repre¬ 
sentations.  The  answer  depends  on  the  type  of  smoothness  considered  and  on  which 
representation  is  used.  The  reader  is  referred  to  Ball  [13]  for  a  thorough  discussion 
of  the  relationship  between  the  smoothness  of  if),  if),  and  if).  We  will  discuss  some  of 
his  results  below. 

The  easiest  case  is  when  “smooth”  is  interpreted  as  continuous.  For  this  case  we 
have  the  following  simple  result.  Contmuity  of  a-ny  o‘ne  of  the  five  functions  in  (5.1) 
implies  continuity  of  the  others.  To  see  this,  note  that  the  arguments  of  any  one  of  the 
functions  if),  if,  if)^  can  be  expressed  as  polynomials  in  the  arguments  of  the  other  two 
functions,  so  continuity  of  one  of  these  functions  implies  continuity  of  the  other  two. 
Similarly,  since  the  principal  invariants,  for  example,  can  be  expressed  as  symmetric 
polynomials  in  01,02,03,  continuity  of  ^  implies  continmty  of  if).  Conversely,  since 
the  unordered  triple  (01,02,03)  of  eigenvalues  of  A  is  a  continuous  function  of  the 
principal  invariants  of  A,^'^  continuity  of  if)  implies  continuity  of  if).  Since  the  principal 
invariants  and  the  moments  of  A  are  continuous  fimctions  of  A,  continuity  of  any 
one  of  the  functions  •f),if),ip^  implies  continuity  of  if).  Finally,  choose  any  orthonormal 
basis  {61,62,63}  for  V.  Then 

^(ui,  cj,  03)  =  ^(0161  (S)  61  -1-  0262  (S)  62  -I-  <1363  <8)  63) ,  (5.2) 

so  continuity  of  if)  implies  continuity  of  if).  Therefore,  Theorems  1.1,  1.3,  and  2. 1-2.3 

^®Cf.  Truesdell  k,  Noll  [4,  §10]  for  a  proof  of  the  representations  The  representation 

follows  from  if  and  the  identities  (2.2)  and  UIx  =  • 

^^As  noted  by  Serrin  [2,  p.  463]  and  Ball  [13,  p.  709],  this  follows  from  the  well-known  fact  that 
the  roots  of  a  polynomial  are  continuous  functions  of  the  coefficients,  the  polynomial  in  this  case 
being  the  characteristic  polynomial  of  A.  The  roots  are  analytic  functions  of  the  coefficients  on  any 
domain  in  which  the  roots  remain  distinct. 
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also  yield  sufficient  conditions  for  the  continuity  of  the  functions  in  the  representations 
of  the  coefficients  a,  7,  and  6.^^ 

The  relationship  between  the  smoothness  of  a  scalar-valued  isotropic  function 
xj}  on  U  and  the  smoothness  of  the  corresponding  functions  in  the  representations 
(5.1)  is  more  complicated  when  “smooth”  denotes  some  degree  of  differentiability. 
We  consider  the  simplest  case  first,  namely,  when  ^(A)  is  expressed  as  a  symmet¬ 
ric  fimction  of  the  eigenvalues  of  A:  V’(A)  =  ■0(01,02503)-  The  domain  U  oi  tp 
is  the  set  of  all  (01,02,03)  €  such  that  01,02,03  are  the  eigenvalues  of  some 
A  E  U.  H  (01,02,03)  €  if,  then  (oi,o„„o„)  €  U  for  every  permutation  Z,m,n  of 
1,2,3,  and  A  =  X)f_iO,e,-  (g>  e,-  eU  for  any  orthonormal  basis  {e,}  for  V.  Then 
A5  =  X)f_i(o,  +  ^,)e,(8)e,-  €  U  for  sufficiently  small  since  W  is  open  and  ||A- A^||  = 
(E^-i  Hence  (oi  +  Si,  02  -f  ^2, 03  +  ^3)  ^  ^  for  sufficiently  small  Si,  which  shows 

that  U  is  an  open  subset  of  R^.  Since  (5.2)  holds  for  any  fixed  orthonormal  basis 
{e,}  for  V,  it  follows  that  smoothness  ofxjx  implies  smoothness  of  xj)  for  any  reason¬ 
able  interpretation  of  the  term  “smooth”.  A  looser  but  more  descriptive  statement 
of  this  result  is  that  if  0(A)  is  a  smooth  function  of  A,  then  0  is  a  smooth  func¬ 
tion  of  the  eigenvalues  of  A.^^  Therefore,  Theorems  1.2,  1.3,  and  f.  1—4.3  also  yield 
sufficient  conditions  for  the  coefficients  a(A),  j3(A),'/(A),  and  ^(A)  to  be  regarded 
as  (r  =  1,2)  functions  of  the  eigenvalues  of  A.  In  the  other  direction.  Ball  [13, 
Theorem  5.5]  proved  that  0  is  C'  then  0  is  C*-,  for  r  =  0, 1,2.  He  conjectured  but 
could  not  prove  that  this  holds  for  any  positive  integer  r.  However,  his  Theorem  5.7 
does  imply  the  weaker  result  that  ifxf  is  C’'+i  then  0  is  C'^,  for  any  positive  integer 
r;  hence  0  is  C°°  if  0  is  C°° . 

The  relationship  between  the  smoothness  of  0  and  0  is  dramatically  different. 

A 

The  domain  of  0  is  the  set 

W  =  {(4,HA.,ff4):A€W}cm.  (5.3) 

Then  U  =  U  W2  L,  where 

w„  =  {(4,J'Ja.J'/4):  Ae«„)cR».  (5.4) 

Similarly,  let  denote  the  set  of  all  points  {0-1,02,  a  f)  E  U  such  that  n  of  the  a,- 
axe  distinct;  equivalently,  (01,02,  <^3)  ^  iff  there  is  an  A  €  with  eigenvalues 
01,02,03.  Clearly,  is  an  open  subset  of  R^.  It  can  be  shown  that  is  also  an 
open  subset  of  R^;  indeed,  Z4  is  the  interior  of  II ,  whereas  ZYj  U  ZY2  is  a  subset  of  the 

^®Serrin  [2]  and  Man  [10]  proceeded  in  the  other  direction.  In  their  proofs  of  the  continuity  of  the 
coefficients  a,/?, 7,  they  regarded  these  coefficients  as  symmetric  functions  of  the  eigenvalues. 

i®It  is  interesting  to  note  (as  did  Ball  [13,  p.  701]  in  the  same  context)  that  the  eigenvalues 
themselves  (either  regarded  individually  when  arranged  in  ascending  order  or  regarded  as  an  un¬ 
ordered  triple)  fail  to  be  differentiable  at  any  A  €  Sym  for  which  <  3,  i.e.,  at  any  point  in 
Synii  U  Sym2;  cf.  Kato  [15,  §11.6.4].  Thus,  if  V"  is  smooth  then  ip  is  the  composition  of  a  smooth 
function  0  :  Z?  C  R®  -+  R  with  the  nondiffereniiable  eigenvalue  function  A  >-*  (01,02, 03). 
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boundary  of  Thus,  unlike  U  and  the  domain  W  of  ^  is  not  an  open  set  (unless 
U  =  i/3).  Since  the  eigenvalues  of  A  axe  analytic  functions  of  the  principal  invariants 
of  A  on  the  open  set  ZY3,  the  same  arguments  used  for  the  case  of  continuity  imply 
that  xb  is  smooth  on  U%  iff  ^  is  smooth  on  iff  xp  is  smooth  on  However,  there 

A  A  A 

can  be  a  loss  in  the  smoothness  of  xp  at  points  in  i/j  U  U^-  One  complication  which 
must  be  addressed  in  discussing  the  smoothness  of  ^  is  the  meaning  to  be  assigned  to 
the  term  “smooth”  at  these  boundary  points.  For  simplicity  we  now  consider  the  case 
when  U  =  Sym-,  then  UZV2  is  the  entire  boundary  ofU.  Ifxj)  is  onU  =  Sym  ( and 
hence  ^  is  C^’’  onU  =  then  xj^  is  C'*'  on  U  in  the  sense  that  there  is  a  C’’  function 
f  on  R3  xvhose  restriction  to  U  equals  In  particular,  if  xj)  is  onU  =  Sym,  then 
^  is  on  U,  i.e.,  ^(A)  may  be  regarded  as  a  continuously  differentiable  function 
of  the  principal  invariants  of  A.  However,  the  requirement  that  ^  be  cannot  be 
weakened  to  For  suppose  that  xfp  has  the  form 

V’(A)  =  xl>{ai,  02, 03)  =  ff(ai)  +  5(02)  +  ^(03) ,  (5.5) 

where  ^  ;  R  — >  R  is  but  not  C^.  Then  xj;  is  C^,  and  hence,  by  the  theorem  of  Ball 
discussed  in  the  previous  paragraph,  xp  is  C^.  But  by  Theorem  3.8  in  Ball  [13],  the 
fact  that  g  is  not  implies  that  for  xp  of  the  form  (5.5),  the  derivative  of  xp  cannot  be 
extended  continuously  from  the  interior  of  U  to  the  boundary  of  ti,  and  hence  xp  fails 
to  be  continuously  differentiable  (in  any  reasonable  sense)  at  the  boundary  of  U.  In 
other  words,  ifxp{A)  is  only  a  function  of  A,  then  xp  may  fail  to  be  a  continuously 
differentiable  function  of  the  principal  invariants  of  A  at  values  of  I^,  Haj  -^a. 
xvhich  A^.  =  0,  i.e.,  for  xvhich  *A  <  3. 

Note  that  none  of  the  theorems  in  the  previous  sections  yield  sufficient  conditions 
for  the  coefficients  a,  /?,  7,  or  0  to  be  C^.  Thus,  none  of  our  smoothness  theorems 
guarantee  that  a,  7,  or  6  can  be  regarded  as  continuously  differentiable  functions  of 
the  principal  invariants  of  A.  Of  course,  since  Proposition  4.1  is  valid  for  arbitrary  r, 
it  might  be  possible  to  apply  the  techniques  of  Section  4  to  obtain  sufficient  conditions 
for  the  continuity  of  higher  order  derivatives  of  the  coefficients.  For  r  =  1,2,3,  we 
have  shown  that  if  $  is  C’’  on  U,  then  6  is  C''~^  on  U,  and  a,  /?,  7  are  C''~^  on  U2  U  U3 
and  on  U.  It  is  reasonable  to  conjecture  that  this  result  also  holds  for  r  >  3. 
However,  a  glance  at  the  formula  (4.2)  for  reveals  that  the  formula  for 

will  be  quite  messy,  and  it  is  not  clear  that  the  latter  formula  can  be  solved  for  D^9 
or  Thus,  even  for  the  next  simplest  case,  r  =  4,  the  approach  developed  here 
might  not  be  useful. 

^°Cf.  Ball  [13,  Lemma  3.2]  for  this  and  other  relations  between  the  topological  properties  oiU  and 
U.  That  Ui  U  U2  lies  in  the  boundary  of  U  also  follows  from  the  fact  that  UIa  must  satisfy 

the  inequality  Aa  >  0  (cf.  (1.13)),  with  Aa  =  0  iff  A  G  //i  UZ/2- 

^^Ball  [13]  attributes  this  result  to  Barbancon  (1972).  The  result  also  follows  from  Theorem  3.2 
in  Ball’s  paper,  which  applies  to  more  general  (but  not  arbitrary)  open  invariant  subsets  U  of  Sym; 
cf.  the  discussion  on  pp.  700  and  705-706. 
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6  Proof  of  Proposition  4.1 

Our  proof  of  Proposition  4.1  employs  Propositions  6.1  and  6.2  below.  We  begin 
with  some  definitions  and  results  which  are  useful  for  a  concise  statement  and  proof 
of  Proposition  6.1.  Let  Orth  denote  the  set  of  orthogonal  tensors.  The  orbit  of  a 
symmetric  tensor  A  is  {QAQ^  ;  Q  €  Orth}.  Thus  the  orbit  of  A  is  the  smallest 
invaricint  subset  of  Sym  which  contains  A;  and  C  Sym  is  invariant  iff  contains 
the  orbit  of  every  A  gU.  By  an  orbit  of  the  invariant  set  U  we  mean  a  subset  of  U 
which  is  the  orbit  of  some  A  gU.  Then  U  is  the  union  of  its  orbits,  and  each  orbit 
of  U  lies  in  one  the  the  invariant  subsets  UxMi-Mz-  A  scalar- valued  or  tensor- valued 
function  on  U  is  isotropic  iff  its  restriction  to  each  orbit  is  isotropic.  In  particular,  a 
scalar- valued  function  on  U  is  isotropic  iff  it  is  constant  on  every  orbit  of  U. 

Consider  a  scalar- valued  isotropic  function  V’  Let  Tr{Sym)  denote  the  space 

of  (covariant)  tensors  of  order  r  on  5t/m,  i.e.,  the  space  of  all  multilinear  maps  from 
Sym  X  •  •  •  X  Sym  (r  times)  into  R.  If  V’  is  r-times  differentiable  at  the  point  A  gU, 
then  D’^rp{A)  G  %{Sym).  Let  U  denote  the  set  of  all  points  in  U  at  which  tp  is  r-times 
differentiable,  and  let  ^  %{Sym).  Then  it  is  not  hard  to  show  that  U 

is  an  invariant  subset  of  U  and 

«^,(QAQ")[QElQ^ . . . ,  QE,Q^]  =  9,(A)[E,, . . . ,  E,]  (6.1) 

for  each  A  gU  and  each  Q  €  Orth;  here  and  below,  Ei,...,Er  denote  arbitrary 
symmetric  tensors.^^ 

Now  consider  any  invariant  subset  IJ  of  Sym  and  an  arbitrary  function  ^ 
%{Sym).  We  say  that  is  isotropic  on  U  if  (6.1)  holds.  This  is  consistent  with  our 
previous  use  of  the  term  for  functions  from  U  into  Sym.  To  see  this,  note  that  any 
B  G  Ti(Sym)  (the  space  of  all  linear  functions  from  Sym  into  R)  can  be  identified 
with  a  imique  B  G  Sym  via  the  canonical  isomorphism  B[E]  =  B*E.  Thus  any 
function  :  U  — >  Sym  may  be  identified  with  a  function  :t{  Ti{Sym).  And  if 
the  latter  function  is  isotropic  in  the  sense  of  (6.1),  we  have 

Q"^i(QAQ")Q.E  =  «^a(QAQ").QEQ"  =  ^i(QAQ")[QEQ^] 

=  ^i(A)[E]  =  ^i(A).E. 

Since  this  holds  for  every  symmetric  tensor  E,  it  follows  that  Q^'9?i(QAQ^)Q  = 
(A)  and  therefore  '®i(QAQ'^)  =  Q^i(A)Q^,  i.e.,  is  an  isotropic  tensor- valued 

function  as  defined  in  the  Introduction.^^  Similarly,  if  we  set  7^(5i/m)  =  R  and 

^^Also,  if  =  D’'ip  then  '®'r(A)  is  a  symmetric  tensor  of  order  r  on  Sym,  i.e.,  $r[Ei, . . . ,  Er]  = 
^r[E<,(i), . . . ,  E,7(r)]  for  any  permutation  a  of  1, . . . ,  r;  cf  (8.12.4)  in  Dieudonne  [9].  However,  this 
property  of  is  not  used  for  any  of  the  results  in  this  section. 

similar  argument  shows  that  if  U  is  the  set  of  all  points  where  the  isotropic  tensor- valued 
function  #  :  W  ^  Sym  is  r  -  1  times  differentiable,  then  may  be  identified  with  an  isotropic 

function  :  Z?  — »  Tr(Sym). 
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interpret  (6.1)  for  r  =  0  as  the  condition  ^o(QA.Q^)  =  '9?o(A),  then  '9q  is  a  sceilar- 
valued  isotropic  function. 

From  (6.1)  it  follows  that  :L(  %.{Sym)  is  isotropic  iff  its  restriction  to  each 

orbit  of  U  is  isotropic.  By  setting  Q  =  and  A  =  RBR^  in  (6.1),  we  see  that 
isotropy  of  implies 

«',(RBRy  )[Ei, . . . ,  E,]  =  ^r(B)[R''  E,R, . . . ,  R?"  E,R]  (6.2) 

for  any  B  €  ZY  and  any  R  €  Orth.  Conversely,  if  for  every  orbit  of  ZY  there  is  at  least 
one  tensor  B  in  the  orbit  for  which  (6.2)  holds  for  each  R  €  Orth,  then  is  isotropic. 
For  if  A  is  any  tensor  in  the  orbit  of  B  then  A  =  QBQ^  for  some  Q  E  Orth,  so  for 
any  Qe  Orth  we  have 

«,(QAQ")(...,QE,Q",...]  =  «,((QQ)B(QQ)")[...,QEiQ',...] 

=  ®,(B)[. . . ,  (QQ)'(QEtQ^)(Qq), . .  .j 

=  9>,(B)|...,QXQ,.-1 

=  ®,(QBQ’')[...,Ei,...] 

where  the  second  equality  follows  from  (6.2)  with  R  =  QQ  and  E*.  — >  QEjtQ^,  and 
the  fourth  equality  follows  from  (6.2)  with  R  =  Q.  Given  the  value  of  at  some 
point  B,  we  can  use  (6.2)  to  extend  fo  an  isotropic  function  on  the  orbit  of  B. 

Recall  that  the  set  ZY<{e,}  C  Z4  H  5{e,}  consists  of  all  A  =  a.e,-  0  e,-  €  Z/  for 

which  Oj  <  02  <  C3. 

Proposition  6.1  If'^r  :  %{Sym)  is  continuous  and  isotropic  on  14^,  then 

the  following  conditions  are  equivalent: 

(1)  ®r  has  a  continuous  extension 

(2)  VB  e  Z/i  U  U2,  the  lim  ’$’r(-A-)  exists; 

A— B 
AeWj 

(3)  VB  €  Z^i  U  U2,  each  principal  basis  {e,}  for  B,  and  any  Ej, . .  .,Er  €  Sym, 
the  lim  'Sr(A)[Ei, . . .  ,E,.]  exists. 

A.--+B 

A€W3nS{e,} 

(4 )  for  each  orbit  of  U  U2  there  is  some  B  in  the  orbit  and  some  principal  basis 
{01,62,63}  for  B  corresponding  to  the  ordered  eigenvalues  bi  <b2  <  63  of  B, 

such  that  the  lim  ^r(A)[Ei, . . . ,  E^]  exists  for  any  Ej, . . . ,  Er  G  Sym. 

A->B 

AQU<{ei} 

When  these  conditions  hold,  unique  and  isotropic  on  H. 


Proof:  As  in  the  proof  of  Proposition  2.1,  the  equivalence  of  (1)  and  (2),  and  the 
uniqueness  ^^5  follows  from  the  fact  that  is  dense  in  U\  and  clearly  (2)  (3) 

(4).  So  assume  (4)  holds.  By  the  linearity  of  (4)  implies  that  lim^-^B  ^r(A), 
A  €  ZY<{e,},  exists.  Let  €  %{Sym)  denote  the  value  of  this  limit.  Extend 

to  an  isotropic  function  on  the  orbit  of  B  by  (6.2)  with  place  of  K  B 

lies  in  the  orbit  of  B  then  B  =  ELi  ®  e,  for  some  orthonormal  basis  {e,}.  Let 
{A„}  (n  =  1,2, . . .)  be  any  sequence  in  converging  to  B.  Let  ELi  «i,nei,n0e,-^„, 
with  Cl  „  <  02, „  <  03, be  a  spectral  decomposition  of  A„;  then  o,-,„  6,-  (e  = 

1,2,3)  by  the  continuity  of  the  ordered  eigenvalues.  Let  Q„  be  the  orthogonal 
tensor  which  maps  ei,„  to  e,-  (i  =  1,2,3).  Then  Q„A„Q]^  =  Ef=iO,,„ei  0  e,-  € 
ZY3<{eJ,  and  since  o.-,„  6,-  it  follows  that  Q„A„Q^  B.  Thus  (4)  implies  that 

^,.(Q„A„Q^)[Ei,...,E,.]  ^,.(B)[Ei,...,E,.]  for  any  Ei,...,E,.,  and  hence  that 

®f,.(Q„A„Q^)  ^r(B).  Now  the  sequence  {Q„}  need  not  converge,  but  since 

Orth  is  compact  there  is  a  subsequence  {Qn^}  such  that  Q„^  €  Orth.  Then 

Q„^A„^Q;;^  ^  QcoBQ^  since  A„,  B.  And  since  Q„A„Q;;  ^  B,  we  also  have 
Q„^A„^Q^^  B.  Therefore  B  =  Q^BQ^,  and 

^,(A„J[Ea,...,E.]  =  ’^r(Qn.A„,Q:j[Q„,EiQ"„^,...,Q„,E,Q:j 

=  ^,(B)([E„...,E,]), 

where  the  first  equality  follows  from  the  isotropy  of  on  ZY3,  and  the  last  equal¬ 
ity  from  the  isotropy  of  011  the  orbit  of  B.  We  have  shown  that  for  every 
sequence  {A„}  in  ZV3  converging  to  B,  there  is  a  subsequence  {A„j^}  such  that 
^,.(A„J[Ei, . . . ,  Er]  converges  to  #r(B)[Ei, . . . ,  E^].  It  follows  that 

lim  ^,(A)[Ei, . . . ,  E,]  =  #,(B)[Ei, . . . ,  E,] 

A— B 
Aews 

for  any  symmetric  tensors  Ej, . . . ,  E^.  Thus  lim^_B  ^r(A),  A  €  ^4,  exists  and  equals 
^r(B).  Since  B  was  an  arbitrary  point  in  an  arbitrary  orbit  of  ZYj  U  Z/2)  condition  (2) 
holds.  Since  is  isotropic  on  each  orbit  of  Hi  U  H2,  it  is  isotropic  on  ZYj  U  i/2-  So  if 
we  set  ^3?  then  is  a  continuous  isotropic  extension  of  to  U.  □ 

As  discussed  above,  is  just  a  scalar-valued  isotropic  function  when  r  =  0,  so 
Proposition  2.1  also  follows  from  Proposition  6.1. Here  we  are  primarily  interested 
in  the  case  r  >  1  with  for  some  scalar-valued  isotropic  function  ^  on 

2"* Note  that  condition  (4)  in  Proposition  6.1  (which  for  r  =  0  is  obtained  by  omitting  the 
El , . . . ,  Er)  is  weaker  than  (4)  of  Proposition  2.1. 
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U^.  Then  Proposition  6.1  yields  sufl&cient  conditions  for  the  existence  of  a  continuous 
extension  of  this  rth-order  derivative  from  ioU.  It  is  tempting  to  conclude  that 
the  existence  of  such  an  extension,  together  with  the  assumption  that  ^  is  on 
U,  implies  that  ip  is  on  U .  This  is  not  true  in  general,  that  is,  it  would  not  be 
true  if  LI3  was  just  an  arbitrary  open  dense  subset  of  The  conclusion  is  true  for 
the  case  considered  here  because  the  set  U  —  D  U2  oi  all  tensors  in  U  with 

at  least  two  equal  eigenvalues  is  sufficiently  nice.  Of  course,  this  requires  some  proof. 
Fortunately,  the  result  we  need  has  been  established  by  Ball.^® 

Proposition  6.2  (BeJl  [13])  Letr  >  1  and  let  ^  be  a  (not  necessarily  isotropic) 
function  from  U  into  R”.  If  ^  is  C’'  on  U^,  md  if  9  has  a  continuous 
extension  to  li,  then  ^  is  C'’  on  II. 

Then  Proposition  4.1  follows  from  Proposition  6.2  (with  ^  ■=  =  1)  and  Propo¬ 

sition  6.1  (with  =  D^ip). 


pointed  out  by  Ball  [13,  §2],  the  Cantor  function  provides  a  counterexample  for  r  =  1. 
^®The  result  stated  here  follows  from  the  much  more  general  Proposition  2.2  in  Ball’s  paper, 
together  with  the  fact  (which  is  proved  on  p.  714  of  his  paper)  that  Sym^  U  Sym2  is  closed  and 
sparse  (cf.  his  Definition  2.1). 
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